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Abstract 

The thermalisation rate for long wavelength fluctuations of the gauge field in the 
abelian Higgs model is calculated from the imaginary part of the self energy. The 
calculation is performed for both the symmetric and symmetry broken phase of the 
theory. 


1 Introduction 

The study of gauge theories at high temperatures has found application in the description 
of a quark-gluon plasma in heavy ion collisions and for the physics of the early universe. In 
recent years, it has become apparent that conventional perturbation theory breaks down 
when considering energy scales much smaller than the temperature. In the case of hot 
QCD, this realisation led to the development of a resummation scheme based on hard 
thermal loops|J^, |2 . 

Recently in a paper by Elmfors, Enqvist and Vilja [JJ, the method of hard thermal 
loops was applied to the electroweak theory in order to determine the thermalisation rate 
of fluctuations in the Higgs held. In this paper, we use the same methods to study the 
thermalisation rate of the gauge held in the abelian Higgs model. This is a hrst step towards 
the calculation of the thermalisation rate for a more realistic theory. The thermalisation 
rate may be used to determine the time at which the gauge bosons associated with a grand 
unified theory will ‘freeze out’ in an expanding universe. This may be of use in the subject 
of baryogenesis in the early universe. 
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In the present paper we calculate the thermalisation rate 7 from the imaginary part of 
the temperature dependent Feynman self energy £ via II 

„ = ImS (^ k ) m 

' 2u y J 

where u is the energy of a given mode. We shall be calculating the thermalisation rate of 
large scale fluctuations (that is at zero 3-momentum). We use the real time formalism of 
thermal held theory^, |j : . [7j since it is straightforward to extract the imaginary part of 
the self energy in this formalism ||. 

The paper is organised as follows. In section 2 we give a brief overview of the abelian 
Higgs model and present the corrections to the propagators and vertices in the hard thermal 
loop approximation. In section 3 we use these results to calculate the thermalisation 
rate both above and below the critical temperature at which the phase transition occurs. 
Finally, we discuss the results in section 4. 

2 The Abelian Higgs model 

The abelian Higgs model consists of a complex scalar held minimally coupled to a 17(1) 
gauge held. It is described by the Lagrangian 

£ = - l -F' w F lw + (£^<F)*(£> M $) - p 2 <T$ - §(<F*$) 2 (2) 

with covariant derivative = d^ — leA M and abelian held strength tensor F^ v = d^A u — 
d u Thh We choose p 2 < 0 so that spontaneous symmetry breaking occurs. Expanding <F 
about its expectation value as 

$(z) = -^2 ( u + ^ x ) + 

we may rewrite Eqn.|2| as £ = £ 0 — H(H^, 0, y) where 

£0 = + 

+\d^(t)d^(t) - \m 2 4> 2 + \d^xd^X ~ Md fl xA 11 (3) 

and I7(H M ,0, y) is a potential term containing only cubic and quartic terms. M = ev and 
m 2 = |Au 2 . We assume A ~ e 2 . 

The gauge is chosen so as to simplify the calculation. That is we choose the gauge so as 
to both simplify the explicit form of the propagators used and also to reduce the number 
of diagrams contributing to hard thermal loops. To this end we choose to use the t’Hooft 
gauge (the renormalisable gauge with £ = 1 ). 

We do not use the covariant gauge as the propagators in this gauge must be resummed 
due to mixing between the photon and the unphysical goldstone mode. As shown in [||, 
this resummation leads to real time propagators that contain derivatives of Dirac delta 
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functions. Although there is in principle no problem with using such propagators, it is 
much simpler to use the renormalisable gauges which do not need this kind of resummation. 

The t’Hooft gauge is then chosen as it enables us to reduce the number of diagrams 
contributing hard thermal loops. This is due to the particularly simple form of the photon 
propagator. The only disadvantage of this choice of gauge is that Fadeev-Popov ghosts 
must be included. 

As a result of this choice of gauge, the following terms must be added to the lagrangian. 

Cg.f. = ~ {d^ + Mxj ~ c^n + M 2 + eMcf^j c (4) 

where c is the ghost held. As a result of this choice of gauge, the Goldstone mode acquires 
a mass M equal to that of the gauge boson. 


2.1 Leading order results 


The leading order corrections in the hard thermal loop approximation are known in the 
case of the abelian higgs model ||TT]. [□], [L2|, [L3], [L4], [L5| . We shall not repeat these calculations 
but merely quote the results. At high temperatures, the vacuum expectation value(VEV) 
decreases from it’s zero temperature value, reaching zero at the critical temperature Tq. 
The temperature dependant VEV vg is given by 


Vp = V 


T 2 (2 e 2 

tv3 + 3 t 


(5) 


The masses of the scalar and ghost sectors of the theory, are affected only by the change 
in the VEV. That is the temperature dependant masses mg and Mg are defined by 


m,- 


= J4 


Mg = evg 


( 6 ) 


However for the gauge boson there is a high temperature correction in addition to that due 
to the temperature dependant VEV. This additional contribution depends on the polari¬ 
sation. That is there are three mass corrections corresponding to the spatially transverse 
and longitudinal polarisations and to the 4-vector longitudinal polarisation. These masses 
are given by [ilQ 1 


Transverse: Ur(ko, k) = 


e 2 T 2 


M 2 g + 


hi _i_ 
k 2 + 


k 2 

k 2 


fco In ( ko+IF 
2k l fc 0 —|k| 


Longitudinal: 


n t (4.k) = + 


( 7 ) 


4-longitudinal: 


n ! = 


M 2 
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There are no hard thermal loop corrections to the vertices |T3], [15| . 

We shall be calculating the thermalisation rate at zero momentum. As such we define 
Mt to be the value of the self energy of the transverse and longitudinal modes in this limit. 


Mf. — U T (k 0 , 0 ) — U L (k 0 , 0 ) — Mj + 


e 2 T 2 

9 


( 8 ) 


Above the critical temperature, the symmetry is restored. There is now no distinction 
between the higgs and goldstone fields. The masses for the particles are now 

ffiggs/Goldstone: m 2 = A( T 2 — T£){ |A + 3e 2 } 


Photon: As in Eqns.7-||| but with Mg = 0 


( 9 ) 


3 Calculation of the thermalisation rate 


To calculate the thermalisation rate we need the imaginary part of the self-energy (T^j. In 
the real time formalism, the imaginary part of the Feynman self energy is most easily 
calculated from the 1-2 component of the self energy via the relation B [1 


0\k o 


ImS = 


2 n(|fc 0 |) 


ImS 


12 


( 10 ) 


where n(X) is the Bose-Einstein distribution. We are working in the symmetric version of 
the real time formalism (that is a = 1 in the notation of Landsman and van Weert J3J). 


3.1 Resummed propagators 

To proceed with the calculation we need the explicit form for the 1-2 component of the 
resummed propagators. For the scalar and ghost sector of the theory, this is easy due to 
the simple form of the hard thermal loop mass correction. However care must be taken 
with the photon propagator due to the complicated momentum dependance of the self 
energies. The 1-2 propagators for the higgs and goldstone modes are given by 

Higgs: iA^ 2 (k) = 2ire 2 n(\k 0 \)5(k 2 — m 2 0 ) (11) 

Goldstone: lA f 2 (k) = 2ne 2 ° n(\k 0 \)5(k 2 — Mj) ( 12 ) 

We must take care with the photon propagator since for < k 2 the self energies of the 
transverse and longitudinal modes acquire an imaginary contribution. By dividing the 
propagator up into it’s three polarisations we can treat each case separately. 

1 C H uv 

iA^(k) = P^iA^k) + P£ v iA \ 2 + -p- %A\ 2 {k ) (13) 
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where P^ u and Pj 11 ' are the transverse and longitudinal projection operators. In the rest 
frame of the heat bath they are given by 


P°°(h) = P°‘(h) = Pfi(k) = 0 

Pk r 




P¥(k) = 
P?(k) = (<T 


k 2 

kJ‘k u 

~P r 


P£ v (k) 


(14) 


From this we find that A f 2 (k) = —Af 2 (k). To calculate Af 2 and Af 2 we use the relation 


Ai 2 = e^nd^o^Af - aA 


(15) 


where in Ap we replace ho by ho + is (ho — *e) when ho is greater than (less than) zero. 
Using this we find that the transverse propagator is given by 


/3| fc n I 


?Af 2 (h 0 ,k) = —27re 2 n(|ho|)c>(h 2 — Ihr(h 0 , k)) 

2 M^c (ho, k) 


= -27re /3l 2° l n(|h 0 |)-- 

(h 2 -Re(n r (h 0 ,k))) +M* c (ko,k) 

where IIt is given by EqnJ^ and Mtc is 

e 2 T 2 | h 0 


h 2 > k 2 


hp < k 2 


M TC (h 0 ,k) = 7T- 


1 - 


k 2 


12 l k , 

Similarly, the longitudinal propagator is given by 
?Af 2 (h 0 ,k) = -2yre £1 5 fll n(|ho|)5(h 2 - n L (h 0 ,k)) 

2M 2 c (h 0 , k) 


= —27re^2^ ro(|h 0 |) 


h 2 -Re(n L (h 0 ,k)) + M 4 G (h 0 ,k) 


h 2 > k 2 


h 2 < k 2 


with Ii L given by Eqn.[7] and M L q is 


Ml c (k 0 , k) = - 7 r 


FT 2 (ho| 

6 |k| 


h 2 

1 _ 

k 2 


(16) 


(17) 


(18) 


(19) 


3.2 The calculation 

There are only two diagrams contributing to X 12 . They are shown in Fig.|I]. From these 
diagrams, we End that Iq 2 is given by 

Ei 2 (p Q , p) = —te 2 J d 4 h|4Af|zA( l 2 (h) + (p — 2h) M (p — 2h)^zAf 2 (h)|?A^(p — h) (20) 
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Figure 1: Feynman diagrams contributing to IrnE^ (Dot denotes hard thermal loop cor¬ 
rected propagator) 


Using this we may now calculate the imaginary part of the transverse and longitudinal 
self-energies. 


Im£f 2 (p 0 ,p) 


^ p i(p) ImS i2(Po,p) 

e 2 r d 4 k 
J (27t) 4 


~ k)^P^{k)tAj 2 (k) + P£ u (k)iA^ 2 (k)^ 

iA* 2 (p-k)iA f 2 (fc). 


( k^k u 'l 

■^(P){(P - 2fcr(p - 2fc)" - 4M 2 ^} (21) 


We have divided up the integral in Eqn.^l] into two sections. The first term contains the 
contribution from the physical modes of the photon propagator. The second term contains 
the contribution from the unphysical modes of the photon propagator and the contribution 
from the Goldstone propagator. Evaluating the second term in Eqn.^1] (which we shall label 
as £f 2 ), we End 


£12 = 2e 2 


d A k 


iA^ 2 (p 


k)iA 


? 2 (k) 


(k 2 ~ Mj) 


k 2 


(P-k) ; 

P 2 


= 0 


( 22 ) 


due to the Dirac delta function in Af 2 (k). We therefore End that the only contribution to 
ImS'T’ comes from the physical modes of the theory. This suggests that our calculation is 
gauge invariant. 


ImEr = 


= e 2 M| 


e 2 


^ r d 4 k 


n(\po\) J (2?r; 


-iA 


4 ~f 2 (P - k)PF(p){ PlikyA^k) + P^{k)iA 


12 (k) 


(23) 


Similarly, when calculating ImE^, the unphysical contributions cancel out. 

We now evaluate Eqn.^3] at p 0 = M T . p = 0. The integration splits up into two regions 
corresponding to the two regions over which A^ 2 and A[ 2 are defined (see Eqns. [T^ - [T9|) . That 
is we split up the integral into the regions k 2 > 0 and k 2 < 0. There is no contribution 
from the region k 2 > 0 since at zero momentum we cannot simultaneously satisfy both 
Dirac delta functions in Eqn.[TT|. 
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The region k 2 < 0 does contribute. Performing all but the final momentum integration 
we find 


IrriEr = 


e 2 M 2 


dk- 


k 2 n[\M T - u\)n(u) ^ mT _^ +Ld _ Mrr) 
37r 2 Jk>\M T -u>\ to ti(Mt) 

2Mtc{Mt — to, k) 

{(M T - to) 2 - k 2 - n t (M t - to, k)j 2 + M* c (M t - to, k) 
(Mt — to) 2 — to, k ) 


(24) 


(M T - to) 2 - k 2 


({M t - to) 2 -k 2 - n l (M t - to, k)) + Mj jC {M T - to, k) 


Similarly ImE^ may be evaluated. At zero momentum we find ImE^ = ImE fj . This is 
to be expected since at zero momentum one cannot distinguish between transverse and 
longitudinal polarisations. 

We must now consider the high temperature limit. Examining the temperature depen¬ 
dant terms from Eqn.f23|, and using the Mellin transformation||. |T7| to extract the high 
temperature behaviour we find 


n(\M T - 

n(-Mj') 


j 1 + n(to) + n(M T — to) M T > to 
| n(to — M t ) — n(to) M t < to 
M t T 

uo\Mt — to | 


We therefore find that to 0(e 2 T), the thermalisation rate is given by 


(25) 

(26) 


7(P = 0) = < 27 > 

where A(A~, ^) is a dimensionless function of the ratio ^ and the ratio of the temperature 
to the critical temperature. We have chosen the factor of 247T in Eqn.|27| so as to compare 
the result with the thermalisation rate above the critical temperature (See section 3.3). 


^,4) = 

T c e 2 


r - 2 M 2 

/ dk—&. 

I k>\M T —0)\ 3c o 2 k 


(28) 


k 2 - (M t - toy 


(M t - to) 2 - k 2 - n t (M t - to, k)) + M* C (M T - to, k) 

(M t - to) 2 


(M t - to) 2 - k 2 - n l (M t - to, k)) + M| C (M T - to, k) 


where the overline denotes that the variable has been scaled by a factor of eT . That is 
k = A;. Note that A(A~, ^) is positive dehnite. 
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This expression depends on the validity of the hard thermal loop approximation. We 
are assuming that all the masses are smaller than the temperature by at least one factor 
of the coupling constant e. That is we should only trust the results for ^ < 1. We should 
also be careful that the constant A(^, 4) > e. In addition to this we must also be sure 
that we do not vary the value of 4 so as to invalidate our assumption that A ~ e 2 . The 
temperature dependance of A (^4, 4) is shown in Fig.^j for a range of values of 4- It can 
be seen that near the critical temperature A(A~, 4) ~ 1- However as the temperature 
decreases, A(A-, 4) becomes smaller. This is due to the increasing size of Mg reducing the 
effect of Landau damping. 

One might expect that since A(A-, 4) oc M|, we should expect that A( 4) vanishes 
at the critical temperature. However, this is not the case. In Fig.[^ we show the contribution 
to A(A~, 4) from the two terms in Eqn.|28|. These correspond to the contributions from the 
transverse and longitudinal parts of the photon propagator. It can be seen that near the 
critical temperature the longitudinal contribution dominates. This is due to the behaviour 
of n T and n L near the light cone. As we approach the critical temperature, becomes 
small around the light cone and can compensate for the smallness of Mg in the integrand 
of Eqn.[28|. However, ng- does not vanish and so the contribution to A(^, 4) from the first 
term in Eqn.[28| vanishes as we approach the critical temperature. 


3.3 Above the critical temperature 


In the symmetric phase, only the second Feynman diagram in Fig.|l] contributes to S 12 . 
From this we find that the imaginary part of the self energy at zero momentum is given by 

r d 4 k 

IuiEtOo, p) = -e 2 2 ^|^|) J J^yi P T U (p)\P - 2 %[p - 2k }j A u{p ~ k )^ A u( k ) (29) 

Evaluating this we find that the thermalisation rate at zero momentum above the critical 
temperature is given by 

e 2 T / 4m 2 \ I 

7(p = 0) = 2s( 1 “ -m) Mr>2m (30) 


where Mt = is the value of the transverse or longitudinal self energy at zero momentum 
and m is defined by Eqn.|[ For Mt > 2 m, this diagram does not contribute to the 
imaginary part of the self energy and the thermalisation rate is of a higher order in the 
perturbation expansion associated with hard thermal loops. Solving the equation Mt = 2 m 
we find that Eqn.[30] is valid in the temperature region 


1 < 


T 


T, 


< 


c 


1 2A + 9e 2 
2A + 8e 2 


(31) 


It should be noted that the value of the thermalisation rate at the critical temperature is 
independent of whether we approach the critical temperature from above or below. The 
decay rate as a function of temperature is shown in Fig.^|. This result is consistent with 
the thermalisation rate for scalar electrodynamics 








4 Conclusions 


We have managed to calculate the thermalisation rate at zero momentum for the gauge 
boson in the abelian Higgs model. The next step is to consider extending this work to a 
more realistic theory that contains fermions and utilises a larger gauge symmetry. While 
studying such a theory is more involved, there are in principle no additional problems in 
calculating the thermalisation rate. 

At the order to which we work a second order phase transition is indicated. However 
near the phase transition vp is small and higher order corrections become important. The 
quantitative results may change but more importantly the qualitative picture may change 
and the theory may have a first order phase transition. However we may trust our results 
away from the critical temperature. 
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